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ABSTRACT 

Even though the dark-matter power spectrum in the absence of biasing predicts a 
number density of halos n{M) oc (i.e. a Schechter a value of —2) at the low- 

mass end [M < IO^^AIq), hydrodynamic simulations have typically produced values 
for stellar systems in good agreement with the observed value a ~ — 1. We explain 
this with a simple physical argument and show that an efficient external gas-heating 
mechanism (such as the UV background included in all hydro codes) will produce a 
critical halo mass below which halos cannot retain their gas and form stars. We test 
this conclusion with GADGET-2-based simulations using various UV backgrounds, 
and for the first time we also investigate the effect of an X-ray background. We show 
that at the present epoch a is depends primarily on the mean gas temperature at 
the star- formation epoch for low-mass systems {z < 3): with no background we find 
a ~ —1.5, with UV only a ~ —1.0, and with UV and X-rays a ~ —0.75. We find the 
critical final halo mass for star formation to be ~ 4 x IO^Mq with a UV background 
and 7 X lO^M© with UV and X-rays. 

Key words: galaxies: dwarf — galaxies: formation — methods: numerical 



1 INTRODUCTION 

The mass spectrum of dark-matter (DM) halos collapsing 
from an initial gaussian-perturbed density field has long 
been an object of study; the forma l ism m ost widely used 
today is that of iPress fc Schechtej l|l974l '). who derived a 
number density of the form 



exponential cutolT at a characteristic luminosity L,. If star 
formation traced the underlying dark matter distribution ex- 
actly (i.e. a constant mass-to- light ratio), we would expect 
to observe a ~ —2 at the faint end. However, galaxy surveys 
such as the SDSS have consistently sho wn a dilferent slope , 
in the range o f -1.2 < a < -0. 9 Csee 'Blanton et a ll 12001 



Table 2 of Al imi fc Courtv02005l . andjPcrcz -Gonzalez et al.l 



n{M)dM = M" 



diner 



dM (1) 



dlnM 

where n{M)dM is the number of halos in the mass range 
M to M -I- dM, p is the mean mass density, 8c is the critical 
(linear) density for collapse, and cr is the rms perturbation. 
If we consider small scales, where a is both S> (5c and ap- 
proximately independent of mass, then the Press-Schechter 
form reduces to 

n{M)dM oc M'^dM (2) 

l|Chiu. Gnedin fc Ostrikeill200ll '). 

The luminosity function of galaxies is com monly 
parametrized in the Schechter form (| Schecht eil 1 1 9 76l ) . 

n{L)dL oc {L/L,)"e'^^^*dL- (3) 

i.e., a power-law of slope a at the low-luminosity end with an 
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1200^ . The extent to which galaxies or dark-matter halos 
fail to follow the clustering distribution of linear pertur- 
bation theory is known generally as "bias", and has been 
an active area of phenomenologi cal and theoretica l research 
(see, e.g., the introdu ction tO iBenson et al.ll2000l and §6 of 
ICoorav fc Shethll2002l for an overview): most relevant here 
is the relative bias between galaxies and DM halos, i.e. 
the extent to which galaxies fail to follow the underlying 
DM. At the low-mass end, the discrepa ncy has come to b e 
called the "missing satellites problem" (|Klvpin et al.|[l999l ). 
since naive Press-Schechter theory predicts ~ 300 satellite 
gal axies for the Local Group, while only ~ 40 are observed 
bimon fc GehallioOTl ) — although this number continues to 
increase as data from e.g. SDSS are analyzed. 

One might expect that this relative bias would be a 
concern for those performing galaxy simulations, and indeed 
DM-only simulations which assign gala xies to DM halos af - 
ter the fact (e.g., the Millennium Run of lSpringel et al.ll2005l ) 
do use an ad-hoc prescription to reduce the baryonic mass in 
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small halos, or overpredict small halos (|De Lucia fc BlaizotI 
I2OO7I ). although agreement is im proving as semia nalytic 
models improve in sophistication ()Guo et al.ll2010bh . How- 
ever, hydrodynamic simulations seem to require no such pre- 
scription: they most often report a faint-end mass spectrum 
slope in the vicinity of -0.9 7 -1-2 (IChiu. Gnedin fc Ostrikeil 
I2OOII : iNagamine et alJlioO^; lAlimi fc Courtvl |2005|) While 



there are some exceptions l|Oppenheimer et al.ll2009l report 
a — —1.45 for their best-fit model: Icra in ct al.' 2009 find 
a — —1.96), it would appear this is a result which is diffi- 
cult to get wrong. 

iDekel fc Silkl i 19861 ) were the first to propose a mecha- 
nism of gas being heate d (by local s upern ovae) and driven 
out of low-mass halos. lEfstathiou I l| 19921 ) introduced the 
idea of the ionizi ng background radiation being the relevant 
mechanism, and iQuinn. Katz. fc Efstathiou I (| 19961 ) used a 
smooth particle hydrodynamic (SPH) simulation to show 
that a photoionizing background does strongly inhibit the 
cooling and infall of gas int o halos of virial mass below 4 x 
lO^MfTi. Pfhoul fc WeinberS (1 19961) did t he sam e with a ID 
Lagrangian code. iNavarro fc Steinmet j (|l997l ) showed that 
the UV background reduced the cooled gas accreted to disk 
galaxi es by half, wit h late-accreting gas preferentially af- 
fected. iGnedin I (|2000l ) used simulations of reionization to de- 
termine that the mass scale at which gas is stripped by pho- 
toionization from the background radiation corresponds to 
the "filtering" length scale, the scale at which baryonic mat- 
ter is smoothed compared to the un derlying dark matter in 
linear perturbation theory; however lOkamoto et al.l (|2008l ) 
using GADGET-2 found a somewhat small er critical mass of 
10^° M(7) (similar valu es were obtained bv lHoeft et al. I I2OO6I 
and[c ram et al.ll2007l ) and, significantly, reported that their 
SPH result was well approximated by a simple prescription 
comparin g the gas temperatu re to the virial tempera ture of 
the halo, isenson et alf l|2002') and 'S omervilld (|2002l ) made 
semianalytic calculations of galaxy evolution and likewise 
found that galaxies fainter than I/« have their star formation 
suppressed or "squelched" by a photoionizing background. 
Meanwhile, the role of supernova feedback in driving gas 
from small halos and thereb y suppr essing star formation 
has continued to be studied jYepes et al. 199 ^; Kay et al] 



l2002l : IScannapieco et ai]l2008l ). iPetkova fc Springell (|2010l ) 
have even performed a simulation showing that the UV out- 
put of local stars can affect the SFR. However, while these 
processes are certainly present and important, our focus here 
is on the ionizing background, which we will show to be both 
necessary and sufficient to reproduce the observed Schechter 
a. 

It should be noted that despite the terminology of "crit- 
ical" or "cutoff" masses, a sharp cutoff in galaxy masses 
is not consistent with the observations: there are low-mass 
(dwarf) galaxies, merely fewer than expecte d. Several mech- 
anism s have been proposed to explain this. IShaviv fc Dekell 
(|2003l ) proposed that some intermediate-mass halos (vdic ~ 
25 — 35 km/s, or M ~ 10^" Mq) have star formation at early 
times, then lose their remaining gas at the time of reion- 
ization and become " red and dead" dwarfs at the present. 
iKravtsov et al.l (|2004l ) suggested that present dwarf satel- 
lites such as those around the Milky Way were originally 
larger but were reduced to their present size by tidal strip- 
ping. Such mechanisms may be needed to explain the very 
low luminosity (> 10^ Lq) Milky Way satellites discovered 



in re cent years (|Willman et al.ll2005l : iBelokurov et al.|[200^ . 
I2OO7I) . 

Two re cent studies are of particular interest. 

ISawala et al. I (|2010l ) studied the evolution of dwarf 
galaxies in isolation and found that both feedback and 
a UV background are necessary to reproduce the ob- 
served properties of the Local Group dwarf spheroidals. 
In particular, they saw that the UV background efficiently 
suppresses star formation after th e epoch of reionizat i on for 
systems below 10^ Af©. Similarly. lOkamoto fc FrenkI (|2009| ) 
simulate a Milky Way-like galajcy and its satellites and find 
a threshold circular speed of 12 km/s ( ~ 8 x 10* M0) for 
halos which can form stars. 

We als o note that even though it has been known at 
least since iMadau fc Efstathioul l| 19991 ) that including the 
X-ray background has a substantial effect on the temptera- 
turc of the IGM (and hence, as we will see, on star forma- 
tion), few sim ulati ons have incorporated this component. 
[Ricotti et al.l (|2008l ) studied its effect on dwarf galaxies and 
found a critical baryon retention mass, but only looked at 
2: > 8. We remedy this deficit by including an X-ray back- 
ground in some of our simulations to study its effect on small 
galaxies from moderate red shift to the present. 

In a previous paper (|Hambrick et al.l |2009| . hereafter 
Paper I), we explored the effect of the ionizing background 
on the stellar and gas properties of large elliptical galaxies. 
Here, we use a similar set of UV and X-ray backgrounds 
to examine their effect on the low end of the galaxy mass 
spectrum, namely the satellites of those ellipticals. Section 
2 presents a simple physical argument for why any efficient 
gas-heating mechanism will yield a flattening of the low- 
mass slope. Section 3 presents the details of the simulations 
we performed to test this argument; section |4] gives the re- 
sults of the same. Section [5] is discussion and conclusions. 



2 PHYSICAL ARGUMENT 



Our argument, w hich dates back to iRees I l| 19861 ) and 
lEfstathiou I (|l992l ). compares the escape velocity of gas in 
a virialized halo to the sound speed of that gas. The escape 
velocity for a halo of mass M and radius r is given by 



2GM 



(4) 



We relate M and r by fixing the density at the standard 
"virialized" value, 

p = 200p 

where p is the mean density of the universe at a given 
epoch (we calculate the result for z = 0, and assume fijv/ = 
0.3, as in our simulations). Then the escape velocity becomes 



2G ( ^pM^ 



1/3 



(5) 



36.0 



M 



1010 A/q 



1/3 



km/s. 



The sound speed, meanwhile, is given for a gas of tem- 
perature T and mean particle mass m by 
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where we assume the adiabatic index 7 = 5/3 for non- 
relativistic gas, and for m we assume primordial abundances 
of H and He. By equating these speeds, we have a formula 
for a critical mass: 

/ J. \3/2 

Merit = 2.19 X 10^ ( YO^J Mq (7) 

At a given gas temperature, halos below this critical mass 
will have gas thermal speeds greater than their escape ve- 
locities. Hence they will not be able to hold on to their gas, 
and therefore be unable to form stars. For a temperature 
of 2 X 10'' K, corresponding to the peak of the H cooling 
curve. Merit = 2 x 10® Mq. Therefore we expect halos with 
less than this total mass to have very little star formation, 
and at late times, very little baryonic mass of any kind, as- 
suming there is a mechanism to heat the gas to 2 x 10* K 
by the epoch where star formation in these halos would be 
significant. Heating above this temperature will further in- 
crease the cutoff mass. Moreover, since larger systems are 
formed by the hierarchical assembly of smaller ones, we ex- 
pect this baryon deficit to creep to larger-mass halos as time 
passes, effectively increasing the critical mass at the present 
by some factor / > 1. 



3 SIMULATIONS 

To verify the physical argument presented above, we ana- 
lyze several high-resolution simulations that are performed 
with various ionizing (UV and X-ray) backgrounds. These 
simulations were performed with a slightly updated ver- 
sion of the GADGET- 2-based code used in iPaper j : refer 
to that work for additional details. In brief, the code is 
based on GADGET-2, and the halos are selected from a 
(50 Mpc)'^ box with cosmological (ACDM) initial condi- 
tions, and resimulated in a (1 Mpc)'' box centered on each 
halo, with the central (0.5 Mpc)^ using hi gh-resolution DM 
and gas/star particles (for more d etails see lNaab et al ] |2007l : 
IJohansson. Naab fc Ostr ikcr 2009). Since we specifically se- 
lect high-density regions containing massive galaxies, our 
results are not necessarily applicable to other environments. 
To ameliorate the effects of cosmic variance, we choose seven 
galaxies (that is, seven initial-conditions (IC) boxes), with 
final virial masses in the range 1 — 2 x lO'^Af©. The sim- 
ulations do not include optical depth effects, in particu- 
lar the self-shielding of dense star-forming regions from the 
UV background. They do (differently from lPaper J ) include 
feedback from supernovae though not AGN, following the 
two-phase ISpringel fc Hernguistl (|2003l 'l model, and a simple 
prescription for metal-line cooling using coo ling rates cal- 
culat ed by Cloudy (v07.02, last described in iFerland et al.l 
1 19981 ). which presumes photoionization and coUisional equi- 
librium for the gas and metal atoms, and assumes 0.1 solar 
metallicity. Star formation is performed at a fixed density 
threshold calculated by the code such that ~ 90% of the gas 
above the threshold will be in cold phase: this threshold is 



Pcrit = 1.6 X 10~^* g cm~^ or nH,crit = 0.07 cm~^. Star par- 
ticles have half the mass of original gas particles; thus each 
gas part icle can ultimately turn into two star particles. As in 
IPaper j . all simulations were performed with initially lOO'' 
each of SPH (i.e. baryon) particles and DM particles, with a 
gravitational softening length of 0.25 kpc for the gas and star 
particles and twice that for the dark matter particles. Gas 
and star particles have masses in the range 4 — 7 x 10^ Mq, 
depending on the size of the box; the assumed cosmology is 
{Qm, fifc/f^M, 0-8, h) = (0.3, 0.7, 0.2, 0.86, 0.65). 

We run these simulations using a variety of ionizing 
backgrounds: the same set which was used in IPaper j . 
with two additions. The first, which we call No UV, has 
as the name suggests no ionizing radiation at all. The 
ne xt three backgrounds ar e Old UV, an updated version 
of iHaardt fc Madaul (119961 ): New UV, which has the same 
spectrum as Old UV but an intensity (as measured by e.g. 
the H photoionization rate) which declines as [1 + z)~^ 
above the pea k at z ~ 2.2; and FG UV, t he background 
constructed in iFaucher-Giguere et al.l l|2009l ). The final two 
backgrounds add an X-ray component, which has a spec- 
trum taken from Sazonov, Ostriker, fc SunvaevI l|2004 ) and 
a normalization from Gilli ot al. (2003)- The New UV-I-X 
model takes the New UV model and adds this X-ray compo- 
nent with the same redshift dependence, while FG UV-I-X 
adds this component to the FG UV bac kground but with 
the redshift dependence of Old UV (i.e. iHaardt fc Madaul 
|199Q) to more closely replicate the quasar background that 
the X-ray component is meant to model. The quasar back- 
ground peaks at z ~ 2 and becomes negligible for ^ > 3 
(|Hopkins. Richards fc Hernauis3l2007l ). 



4 SIMULATION RESULTS 
4.1 Results at 2 = 

Having performed the simulations, we identify virialized sys- 
tems in both the star and total (stars plus gas plus DM) 
particle populations using AHF, the Amiga Halo Finder 
(KnoUma nn fc Knebe 2009 ). AHF uses a density grid-tree 
method followed by the removal of unbound particles to find 
halos and subhalos, and returns the virial mass of each. The 
virial mass is defined as the total mass of all (relevant) par- 
ticles inside the virial radius, which in turn is defined as the 
radius where the calculated density profile exceeds a thresh- 
old multiple of the critical density determined by cosmology 
and redshift (~ 200); for subhalos the virial radius is instead 
the largest radius of any particle which is gravitationally 
bound to the halo. See the AHF documentation for details. 

Figure [T] shows the mass distribution of stellar systems 
(identified using star particles only). AHF can identify sys- 
tems as small as 1.5 x lO^M©, or about 40 star particles, 
but we adopt 10* Mq as a lower mass limit to ensure com- 
pleteness. The mean number of systems (averaged over the 
seven ICs) in the lowest mass decade (10*-10^Mg) declines 
from 131 with No UV to ~ 43 with the three UV models 
and ~ 20 for the two UV plus X-ray cases, differences of 
36cr and 14cr respectively, while the number in the two high- 
est mass bins (> 10^" Mq) remains nearly constant. As an 
aside, we expect these highest-mass bins will be high com- 
pared to the global spectrum because we specifically select 
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Table 1. Schechter a values and maximum "barren" 
with various backgrounds. 



halo masses 



108 



109 



lO'o 
M. (MJ 



10" 



lO'a 



Figure 1. The average number of stellar systems (galaxies) per 
decade in stellar mass identified by AHF with mass in given range, 
at 2 = 0. The error bars represent Icr statistical errors among the 
7 initial conditions. Also included is an power law for com- 

parison. The number of < 10^ Mq galaxies decreases markedly 
as one adds first UV and then X-ray radiation. 



boxes which contain a massive galaxy; furthermore, systems 
of this virial mass (Mnaio ~ 2 x 10^^ M©) have stellar masses 
whic h are very vulnerable to suppression by AGN feedback 
fe.g. iMcCarthv et al.ir2010h . not included here. 

We also see that the UV-only models (Old UV, FG UV 
and New UV) are essentially identical to each other in mass 
spectrum. The two models with X-rays (New UV+X and 
FG UV-I-X) are also identical to each other, though of course 
quite different from the models without X-rays. This sug- 
gests that the details of the ionizing background above z — 3, 
including the epoch of reionization, are relatively unimpor- 
tant in determining the mass spectrum. 

Using the Bayesian inference method of 
iKellv. Fan, fc VestergaardI (|2008l ). we estimate the Schechter 
parameters a and M*; the results for a are presented in 
Table [1] This method has the advantage of not requiring 
any binning of the data. Values for M, were roughly 
constant for the various background models in the range 
2^09.5-10^^. g^jgQ note that we use stellar mass as a proxy 
for luminosity without considering, e.g., age. Consistent 
with the appearance of Figure [1] a increases dramatically 
from —1.5 for the no-background case to roughly —1 for 
the UV-only models and —0.75 for the models with X-rays. 
Thus we see that the addition of a simple UV ionizing 
background was the only change necessary (with feedback; 
see below) to move the low end of the stellar mass spectrum 
into consistency with observations, while an additional 
X-ray background (which most simulations historically have 
not used) seems to make the slope too flat. 

To directly compare these results with the mechanism 
presented in Sj2l we estimate f Merit (the effective star- 
formation cutoff mass) by the following method. We take 



Name 


a 




log f Merit 








No UV 


-1.52 


± 


0.06 


8.19 + 0.06 


0.47 


+ 


0.02 


Old UV 


-1.01 


± 


0.04 


8.55 + 0.08 


0.37 


+ 


0.02 


FG UV 


-1.03 


± 


0.06 


8.50 + 0.08 


0.37 


+ 


0.02 


New UV 


-1.08 


± 


0.05 


8.63 + 0.05 


0.37 


+ 


0.02 


FG UV+X 


-0.73 


± 


0.04 


8.95 + 0.07 


0.32 


+ 


0.02 


New UV+X 


-0.78 


± 


0.08 


8.78 + 0.04 


0.32 


+ 


0.02 


No UV-F 


-1.64 


± 


0.07 


8.06 + 0.06 


0.48 


+ 


0.04 


FG UV-F 


-1.29 


± 


0.06 


8.22 + 0.07 


0.38 


+ 


0.04 


FG UV+X-F 


-1.16 


± 


0.11 


8.22 + 0.09 


0.33 


+ 


0.04 



The Schechter-o values for the star particles, maximum "barren" 
halo masses, and overall baryon-conversion efficiency with various 
backgrounds. Adding UV and X-ray heating steadily increases a: 
models with UV only are consistent with observations. Similarly, 
the critical mass for baryon retention increases, and the baryon- 
conversion efficiency decreases, when the UV background and X- 
rays are added. The "-F" rows represent simulations with no SN 
feedback; see text for details. 



the mass spectrum of halos identified by AHF from the to- 
tal population of particles (stars, gas and DM), and find 
for each the total stellar mass in the halo. Then starting 
at the low- mass end of the spectrum we consider an 11- 
point moving average of the ratio of total stellar mass (M*) 
to total overall virial mass (Mnaio) among the halos. (The 
number 11 was chosen to balance between mass resolution 
and the inherent scatter in stellar mass to total mass ra- 
tio among small halos, but even using a 41-point average 
does not significantly change the results.) We move upward 
through the mass spectrum until this ratio exceeds some 
fixed value, which we choose to be half of the globa l aver - 
age (^0.16), which i s the s ame definition as iGnedin 1 (|2000l ) 
and lOkamoto et al.l (|2008h . The particular choice of criti- 
cal mass is not very important, as the baryon-conversion 
efficiency is a steep function of halo mass around the criti- 
cal point (see Figure 2 of ,Okamoto et a l. 2008) — even if we 
choose a threshold of one-tenth the global average, f Merit 
falls by only ~ 0.2 dex. Thus f Merit represents the largest 
total mass at which systems have severely restricted star 
formation. Mathematically: 



Y^M^i) 

f Merit = Mnaio (iV) 9 

^ 2^ Mnaio W 



Af-l-5 

N+5 

E Mnaio (j) 

i=N-5 



(8) 

The results are presented in the third column of Table [1] 
We find the critical mass in the No UV case is 1.5 x 10* Mq, 
but when a UV background is added this increases by a fac- 
tor of more than 2 to ~3.6 x 10* M©, and with X-rays by 
anot her factor of 2 to ^ 7.3 x 10* M©, in good agreement 
with l Okamoto fc FrenkI (12009'), who also looked at satellite 
galaxies. They used a background with UV but not X-rays, 
but they did include supernova winds, which should act in a 
similar way as X-rays to blow out gas from small halos; thus 
we would expect our UV+X critical mass to be the most sim- 
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ilar to theirs, which is what we find. lOkamoto et al.l l|2008h . 
meanwhile, report a significantly higher value of ~ 10^° 
using only isolated halos; see the discussion in the next para- 
graph. We find (see §4.2) that No UV has a gas temperature 
at the star-forming epoch of 2 x 10^ K, so comparing that 
cutoff mass to the value of 2 x 10* Mq derived in we 
find that the dimensionless hierarchical-structure correction 
factor / « 0.75, actually less than unit y. This seems to sup- 
port the idea of iKravtsov et al] (|2004l ) that small satellites 
at late times were originally larger, and indicates that sys- 
tems above the critical mass are able to accrete sufficient 
gas to overcome the deficit of accre ted stars . 

Finally, we revisit the finding in |Paper~3 that overall star 
formation is suppressed with the ionizing background. To 
compare across the various initial conditions, we calculate 
the baryon-conversion efficiency, e*, defined as 

_ 

where M* and Mt are the total stellar and baryonic (gas plus 
stars) masses, respectively, in the high-resolution simulation 
box. The results are in the last column of Table [T] These 
efficiencies are higher by a factor of roughly fou r than the 
observed values, recently measur ed bv lGuo et all (j2010al ) to 
be ~ 20% (or ~ 4% as a fraction of total mass rather than 
baryons) for halos of mass ~ 6 x 10^^ Mq, and less for other 
halo masses. Similar results were derived by iMoster et al.l 
l|2010t ) (who found that effi c iency is linear in mass at low 
masses) and lTruiillo-Gome3 (|2010D . 

Figure [2] shows the baryon-conversion efficiency as a 
function of halo mass, using a different definition: 

Afb, implied nfcMHalo Mualo 

where M-^^c is the stellar mass with i n 0.1 virial radii of 
the halo center, similar to lSales et al.l (|20l d'). and the "im- 
plied" baryonic mass Mi,,inipiied = MHE,io{^b/^M) is simply 
the mass of baryons associated with the halo well before 
reioniza tion. Figure [2] is analogous to Figure 5 of IGuo et al.l 
()2010al '). although the results are not directly comparable 
(see below). There is a factor of 1.5 — 5 reduction in star 
formation with the addition of the ionizing backgrounds, 
and that the effect is stronger in lower mass halos. Ifow- 
ever our results are still too high compared to the SDSS 
results by a factor of ~ 2 in the highest-mass bin, and 
more in the lower ones. We ascribe this discrepancy first to 
fact that we are not looking at any field galaxies, but only 
galaxies in a region specifically selected to have an over- 
density and a massive central galaxy. This is clearly indi- 
cated in the figure, since SF efficiency actually increases 
with decreasing halo mass from 1O"-^M0 to IO^^Mq: a 
satellite galaxy will naturally have a smaller virial radius 
than a field galaxy, and hence a smaller DM mass will 
be found by the halo finder, not to mention the enhanced 
star formation from the denser environment ([Okamoto et al. 
|2008| ). The lack of SN- driven winds in our simulations, 
which would otherwise eject gas from low-mass galaxies, 
also makes a significant contribution; the relation between 
mass- loss in win ds and the SFR has been verified obser- 
vation allv tRupke. Veilleux fc SaiidersI [20051 ) and theoreti- 
cally l|Oppenheimer fc Davj|2005i r We also note that the 
other r ecent simulations s hown i n Figure [2l-40kamoto et al.l 
(|2005l ). iGovernato et all l|2007l ). IScannapieco et al.l l|2009r 
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Figure 2. The baryon-conversion efficiency as function of (to- 
tal, viral) halo mass. The scatter points are previous simulations. 
Our results, the connected error bars, are the median baryon- 
conversion efficiencies for 1 decade wide bins in halo mass across 
all 7 ICs; the error bars are the inter-quartile ranges, and there is 
a small horizontal offset for clarity. The solid lines are the models 
with SN feedback and the dashed lines are the models without 
it (which use 3 ICs rather than 7). The addition of the ionizing 
backgrounds reduces star formation by a factor of < 5, depending 
on the mass bin. 



and IPiontek fc Steinmet j (|2009l ) — seem to have too much 
star formation by factors of a few as well, and furthermore 
we are primarily interested in the relative differences be- 
tween the various background models: thus an overall excess 
of star formation can be safely disregarded. As with a and 
f Merit, there is no significant difference between the three 
UV-only models, or between the two UV plus X-ray models. 

To explore the effect of the SN feedback on our results, 
we ran three sets of simulations (the A, C, and E ICs) using 
the No UV, FG UV and FG UV+X backgrounds with no 
feedback. The primary results are given in the last three rows 
of Table [T] The Schechter a values are somewhat steeper in 
each case compared to the runs with feedback, which we 
expect since feedback gives a direct energy injection at the 
site of star formation, but the differential effect of the back- 
ground is much the same, although not quite as strong. In 
particular, the FG UV-I-X-F model produces an a within 
the observed range, showing that SN feedback, while cer- 
tainly present and important, is not necessary to produce 
the correct spectrum of low-mass galaxies. A similar result 
is obtained for /Merit, which increases by roughly 0.2 dex 
from no UV to FG UV and changes negligibly when X-rays 
are added. It therefore appears that the effects of the ioniz- 
ing background and of SN feedback amplify o ne another in a 
sup erlinear fashion, as s ugg ested bv lPawlik fc Schave (20oi) 
and ISawala et all (|201G| ). This is also indicated by Figure 2, 
where the dashed lines correspond to the no-feedback mod- 
els: in the lowest-mass bin, adding feedback causes a greater 
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relative reduction in baryon-conversion efficiency when an 
ionizing background is present. 

T o make a detailed comparison with ISawala et alJ 
we calculate the baryon-conversions efficiencies for 
a 0.2-dex-wide bin around their fiducial halo mass of 7 x 
10* Mq. For their ru ns with UV backgro und only, SN feed- 
back only, and both, ISawala et al.l (|2010l ) give efficiencies of 
0.77, 0.05, and 0.03, respectively. Here, however, the respec- 
tive models FG UV-F, No UV and FG UV giv e median val- 
ues of , respectively, 0.16, 0.47, and 0.04. So for lSawala et al.l 
l|2010h . the feedback is by far the most important contrib- 
utor, whereas here the background is somewhat more im- 
portant. The stronger effect of the bac kground in th i s work 
seems likely to be the effect of timing: ISawala et al.l (|201G| ) 
turn their background on at 2 = 6, by which time their 
galaxy has already created 20 — 30% of its final stellar mass; 
in our simulations, however, galaxies of this final size have 
no significant star formation until z ^ ?> (see the next sub- 
section). Therefore our galaxies show no significant change 
in baryon conversion between Old UV, which turns on at 
z = Q, and FG UV, which turns on aX z = 10, but we specu- 
late that if the lSawala et al.l (|20ld ) runs were to be repeated 
with FG UV, the background would have a more significant 
effect. The relatively weak effect of feedback for our models, 
on the other hand, is likely the effect of environment: the 
relatively dense environment and corresponding IGM pres- 
sure on our small halos means that they do not have enough 
supernova energy to unbind their gas. 

The physical mechanisms governing the interaction be- 
tween feedback and the background radiation are beyond 
the scope of this paper, but our results are consistent with a 
picture where feedback moves gas from the centers of halos 
to the outer regions, where the background then provides 
enough energy to unbind the gas altogether. Thus there is 
no change in the /Merit when X-rays are added in the ab- 
sence of feedback because the dense central gas has too short 
a cooling time to be affected (although the smallest halos 
with the lowest central densities can still be suppressed, in- 
creasing a), whereas if it were first pushed out somewhat by 
feedback the density would be lower, the cooling time longer 
and X-rays could have an effect. Conversely, with feedback 
and no background (the No UV case), the gas gets pushed 
outward, but not enough to unbind, so it falls back in and 
forms more stars. 



4.2 Evolution with Time 

To investigate the origins of the discrepancies in low-mass 
systems, we examine the simulations at higher redshifts. Fig- 
ure [3] shows the stellar mass spectrum of the FG UV simu- 
lations from 2 = 4 to the present. The number of low-mass 
systems peaks at 2; ~ 2 — 3 and then declines as small sys- 
tems merge into larger ones. This peak is som ewhat la ter 
than the overall star-formation peak, found in IPaper j to 
be at « « 4, which is pre c isely the " downsizing" effect first 
reported bv ICowie etall l|l996h fsee lZheng et al.ll2007l . for 
more recent observational results). This effect explains why 
our three UV models show such similar results here: ai z — 4 
the Old UV and FG UV models have not finished Hell reion- 
ization, and Old UV has ~ 3 times lower ionization rate than 
the other two models, while for z ^ 3 the three models are 
nearly identical in ionization rate and have all fully reion- 
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Figure 3. The mass spectrum (average number of systems per 
decade of mass) of stellar systems with the FG UV background at 
various epochs. The number of small galaxies peaks at 2 < z < 3. 



ized (see Figure 1 of IPaperT) . Thus at the epoch when the 
smallest halos are forming stars, the gas temperatures in 
the three models are much more similar than when stars in 
more dens e regions formed (gas temperatures are plotted in 
Figure 2 of lPaper J) . On the other hand, the models with X- 
rays have higher gas temperatures than the UV-only models 
for all 2 < 4, which is refiected in the increased suppression 
of small systems. 

Table [2] compares the total number of identified (stel- 
lar) systems in the whole box for the No UV and FG UV 
cases, as well as the baryon conversion efficiency e^^aH, de- 
fined in Equation 9, for two bins in halo mass (10*'^ < 
(A/Haio/Me) < 10»-5 and lO" '^ < (Mnaio/Mo)). Although 
No UV has ~ 3 times the number of systems as FG UV, 
they both show the same time dependence: the number of 
systems reaches a peak at 2 < 2 < 3, followed by mergers re- 
ducing the number. The stellar mass fraction for high-mass 
halos is nearly the same for the two models, while for low- 
mass halos the difference is a factor of 2 — 6. Interestingly, in 
the FGUV case the low-mass stellar mass fraction declines 
somewhat from 2 = 2 to 0, presumably because closer satel- 
lites of the central galaxy, which had been able to form more 
stars than the average because of their denser environment, 
are finally accreted and lost (while with No UV all satellites 
can form many stars). 

Next, we test the validity of our physical argument by 
comparing the mean gas temperature of the simulations at 
2 = 3 (when the small systems are forming) to their critical 
cutoff masses at 2 = 0. (In calculating the mean gas temper- 
ature we do not separate dense galactic gas from the IGM; 
however at 2 = 3 only ~ 5% of gas is at virial densities, 
so any contamination is less than the variance of our ICs.) 
Our theory predicts Merit oc T^'^. Figure U shows the data 
from our six background models together with the best-fit 
power law. We find that Aferit oc yi-54±o.i7^ good agree- 
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Table 2. Stellar system statistics over time 





low- 


mass 


high- 


mass 






z 


No UV 


FG UV 


No UV 


FG UV 


No UV 


FG UV 


4 


.111 


.052 


.546 


.569 


153 


63 


3 


.266 


.140 


.610 


.618 


425 


159 


2.5 


.325 


.180 


.640 


.650 


503 


174 


2 


.400 


.220 


.646 


.637 


567 


169 


1 


.467 


.165 


.728 


.730 


526 


134 


0.5 


.475 


.118 


.739 


.792 


449 


109 


0.0 


.475 


.080 


.787 


.776 


328 


82 


2.5 


.173 


.353 


.697 


.712 


259 


443 


0.5 




.010 




.804 




166 



The baryon conversion efficiency for low-mass (lO*'^ < 
(MHaio/Mo) < 109-5) and high-mass (IQiO'S < (Mnaio/Mo)) 
halos, and total number of stellar systems, for the No UV and FG 
UV backgrounds at various epochs. Statistical errors are ^ 20% 
for the first six rows. The last two rows, below the line, are the 
200^ simulations of halo A. 




10'' 105 
T (K) 

Figure 4. The critical star-free halo mass at 2: = versus the 
mean gas temperature at 2 = 3 for each background model, to- 
gether with the best-fit power law. Vertical and horizontal error 
bars represent la statistical errors among the 7 ICs. We find 
log Merit <^ (1.54 ±0.17) logT, in good agreement with the theo- 
retical value of 1.5. 

ment with the theoretical value. The goodness-of-fit statistic 
X^/d.i. = 1.8 for this fit. Notice that even in the No UV case 
stellar feedback and dynamical heating are sufficient for the 
gas to reach 2 x 10* K hy z — 3 (the gas has reionized by 
this epoch), but the UV heating increases the temperature 
by a factor of 1.6, and the X-rays by another factor of 1.5. 

Finally, we examine the level of numerical convergence 
in our results. We perform simulations of galaxy A with 
200^ SPH particles (i.e. twice the spatial resolution and 



eight times the mass resolution) with the No UV and FG 
UV backgrounds (and including feedback). Due to compu- 
tational constraints, the No UV simulation was run only 
to z = 2.5, and the FG UV to z = 0.5. At 2 = 2.5, the 
Schechter a values for No UV and FG UV with 200^ reso- 
lution are —1.39 and —1.33, respectively. These values are 
more similar than the 100"^ runs, which have a = —1.45 for 
No UV and a = -1.07 for FG UV at the same redshift. The 
steepening slope for FG UV as resolution increases is ex- 
plained by smaller stellar systems which are allowed by the 
higher peak densities associated with increased mass reso- 
lution, and which merge away at late times to flatten the 
low-mass slope once again (whereas No UV can keep mak- 
ing new small systems). Indeed, at z — 0.5, FG UV 200'' 
has a — —1.10, compared to the 100^ run a = —0.95 (and 
No UV 100^ a = -1.57). Meanwhile, the No UV 200=* sim- 
ulation at z = 2.5 still has a large number of small stellar 
systems which have not yet formed stars, and whose gas is 
disrupted by the background in the FG UV run. In fact, 
65% of the identified halos in No UV 200^ at z = 2.5 are in 
this category, compared to only 39% for the 100'^ case. How- 
ever, the overall baryon-conversion efficiency is not signifi- 
cantly different as a function of resolution, at either z = 2.5 
(< 8% difference) ot z = 0.5 (3%). The last two rows of Ta- 
ble [2] show the breakdown by mass. FG UV 200^ has a much 
lower efficiency than 100^ in the low mass bin aX z — 0.5. 
This is the result of the same mechanism just discussed: the 
increased resolution allows small halos to form and merge 
while remaining star-free because of the UV background. 
Further, at z = 2.5 FG UV 200^ has a higher efficiency 
than No UV 200^. This is again the effect of proximity: 
while No UV has twice the number of halos in the low- 
mass bin as FG UV, FG UV's are preferentially close to 
the central galaxy and thus have lower virial masses and 
more efficient star formation, and increased numerical res- 
olution allows them to exist even closer without being lost 
to the halo finder. (At the 25"^ percentile of galaxies, FG 
UV has roughly a tenth the efficiency of No UV in this mass 
bin at both resolutions.) Finally, the /Merit statistic is not 
strongly affected by the increase in resolution: at z — 2.5, 
log /Merit = (8.54,8.76) for (No UV, FG UV), compared to 
(8.47, 8.85) for the respective lOO"' runs. Thus we assert that 
our results at z = would not be qualitatively different with 
increased resolution, and quantitatively different by at most 
15%. 



5 DISCUSSION & CONCLUSIONS 

In lPaper^ we examined the detailed effects on galaxy prop- 
erties of modest changes to the ionizing radiation back- 
ground (specifically the intensity at high 2). Here we have 
simulated seven different galaxies with their satellites, using 
six different models of the ionizing background, and found 
that any UV background is sufficient to replicate the obser- 
vational results of the low end of the galaxy mass spectrum, 
as parametrized by Schechter a. The details of the spec- 
trum and redshift dependence of the background seem to 
have no significant effect. We find that the number of small 
(< 10^ M© in stars) systems at z = correlates with the crit- 
ical star-free halo mass, which in turn correlates well with 
the mean gas temperature at z = 3. Further, the value of 
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our critical mass agrees with other recent studies of small 
gala^xies. These simulation results are consistent with our 
theoretical picture, where only halos with an escape veloc- 
ity greater than the sound speed of the gas can retain their 
gas long enough for it to cool and form stars. 

Thus we find the addition of an X-ray background 
causes additional suppression of low-mass halos (by further 
increasing the mean temperature of the gas), making the 
low-mass spectrum too flat compared with observations. It is 
possible that this is the result of optical depth/self-shielding 
effects not being included in the code, thus making the UV 
background more efficient at heating/ionization than it is in 
reality. Optical depth in general decreases with harder radi- 
ation, so the optically-thin assumption is valid for X-rays. 
However, the small halos which we study here have the low- 
est densities and t hus are the least effective self-shielders: 
ISawala et all (|2010l ) found that systems below the threshold 
mass of 10^ Mq had their gas densities sufficiently reduced 
by feedback and that self-shielding of the UV background 
was not important. Our resolution study weakly indicates 
that a may decrease with increasing resolution, and a future 
paper (Hambrick et al., 2010, in prep.) finds that includ- 
ing AGN thermal feedback causes a to decrease modestly 
as well, presumably by driving gas back out of the central 
galaxy to form stars in satellites. This is certainly an area 
deserving further study. 

We also study the baryon-conversion efficiencies of small 
halos, and find that for small halos in the neighborhood 
of massive ones, and find that the background is at least 
equally important as SN feedback for lowering the efficiency 
to observed values, in contrast to isolated dwarfs where feed- 
back is the dominant mechanism. 

In sum, reproducing the observed numbers of low-mass 
galaxies in hydrodynamic simulations is no problem at all, 
being adequately resolved by the heating mechanisms al- 
ready included in standard hydro codes: feedback and the 
ionizing background. 

JPO was supported by NSF grant AST 07-07505 and 
NASA grant NNX08AH31G. DCH thanks Ludwig Oser for 
invaluable debugging assistance. TN and PHJ acknowledge 
support by the DFG cluster of excellence 'Origin and Struc- 
ture of the Universe'. 
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